基于光学测角数据的风云四号同步轨道卫星精密定轨

针对风云四号同步卫星的精密定轨和精度评估需求，首先利用地面光学测角数据对FY-4A卫星进行精密定轨，定轨后方位角和高度角的残差rms分别为0.25＂和0.45＂。与基于测距数据的轨道相比，位置精度在有测角数据的弧段内小于50m。进一步联合测角数据和测距数据对FY-4A卫星进行联合定轨，定轨后轨道重叠精度优于15m。利用联合定轨结果评估了基于测距数据的实时轨道产品精度，可以明显发现轨道精度随着测距数据的积累而逐步提高。     

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

椭圆轨道卫星对地全球覆盖电推进控制

针对椭圆轨道卫星近/远地点的星下点对全球或特定纬度区域的访问问题，提出一种连续小推力下的对地覆盖控制策略。首先，推导了自然摄动对卫星拱线变化的影响，并探讨了进行小推力覆盖控制的必要性。然后，针对燃料消耗的优化问题，将控制方程展开成含傅里叶级数的形式，用以获得便于星上计算的解析形式的次优解，同时探讨了截取阶数与优化程度的关系。在进行拱线控制的同时，通过合理设置约束，对椭圆轨道的近地点高度进行保护，确保卫星安全运行。仿真结果表明，提出的方法能够以适当的燃料消耗代价实现椭圆轨道的近/远地点的全球覆盖控制或特定纬度区域的反复推扫，且控制力在可接受的范围内。     

低排放燃烧室NO_x排放试验及预测方法研究

为了研究中心分级贫油低排放燃烧室的排放特性和排放预测方法,针对一个低排放头部方案,在单头部燃烧室试验件上,在不同的温度、压力、油气比、供油模式和分级比条件下,测量其排放性能。以Lefebvre排放经验预测公式为基础,采用经验分析方法拟合排放试验数据,归纳出适用于本头部方案的排放预测公式。表征预测好坏的判定系数R2在小工况下和大工况下分别为0.95和0.93,表明预测结果与试验结果符合度较好。小工况和大工况排放特性不同,对仅预燃级喷油的小工况工作模式,NO_x排放主要受化学恰当燃烧温度和预燃级局部当量比的影响;对预燃级和主燃级同时喷油的大工况工作模式,NO_x排放主要受燃烧区温度和主燃级燃油比例的影响。     

火星轨道交会自主导航与制导方法

对火星采样返回任务中的火星轨道交会自主导航和制导技术进行了研究。采用光学自主导航敏感器测量的火星中心方向和视半径，相对敏感器测量的相对位置等观测量，设计了导航滤波器同时估计轨返组合体和上升器的轨道。在导航滤波器设计中，针对光学自主导航敏感器更新频率远低于滤波解算频率的问题，设计了一种连续观测量构造算法，确保每个滤波周期均可进行测量更新，以提高导航精度。基于导航滤波器估计结果，采用T-H制导设计了4脉冲共椭圆交会策略实施轨道控制，从而构成近程交会自主导航和制导方案用于完成火星轨道交会任务。通过数学仿真校验了所提出方法的有效性。     

基于动态补偿的多电机控制算法

多电机协调控制达到转速一致是电机控制中的一个关键难题，外加负载变动和电机参数变化时会引起电机性能的下降，达不到良好的控制效果。为使电机间保持转速同步，建立了永磁同步电机d-q坐标系下的矢量控制数学模型，基于转速跟随控制，提出一种基于单神经元PID的变增益速度补偿器进行偏差耦合控制。在MATLAB/Simulink建立了3台永磁同步电机的仿真模型。仿真结果表明，与传统PID固定增益速度补偿器算法相比，单神经元PID的变增益速度补偿器具有更强的鲁棒性以及更快的收敛性。     

同步双小行星系统共振轨道设计

研究了同步双小行星系统中共振轨道的设计方法及演化规律。首先，基于双椭球模型建立探测器运动方程，并给出共振轨道初值选取方法。然后，利用改进并行打靶法，提出一种双小行星系统平面共振轨道两步修正方法。同时结合稳定性理论及分岔理论，给出双小行星系统三维共振轨道生成和延拓方法；最后，以双小行星系统1999KW4为例，设计了共振比为1∶1,1∶2,1∶3,1∶4,2∶3的平面和空间共振轨道族，并分析了共振轨道的特性及轨道周期和轨道能量的变化规律。给出的双小行星系统中共振轨道的设计方法具有普适性，对未来双小行星系统探测任务中的轨道设计具有一定的参考意义与借鉴价值。     

部署retro-GEO巡视器的月球借力飞行轨道分析

retro-GEO是指逆行（retrograde）地球静止轨道（geostationary Earth orbit, GEO），该轨道与GEO轨道高度相同或相近，但倾角为180°，安装在retro-GEO卫星上的巡视器可每12h对GEO带空间资产附近碎片环境安全预警。直接西向发射retro-GEO卫星存在地面测控和发射能耗较大的困难。基于平面四体模型，为降低设计变量敏感性，以近月点参数为设计变量，建立了部署retro-GEO巡视器的月球借力飞行轨道设计模型，利用轨道动力学模型延拓策略，得到该类轨道绕月后返回地球飞行时长只能约为114.79h，该结论可用于求解该类轨道高精度轨道动力学模型解。     

轨控大干扰下的系统角动量管理

摘要： 针对利用角动量交换装置系统实现轨控期间姿态维持的情况，研究当轨控推力器存在大力矩扰动时系统角动量管理方法.结合对地运行航天器轨道运动特性，分析轨控干扰力矩的积累角动量变化规律；考虑角动量管理装置系统的角动量存贮容量，基于角动量积累规律提出了一种结合偏置角动量建立与对称分布轨控位置选择的系统角动量管理方法，实现了对称轨控位置积累角动量互相抵消的自平衡效果，极大程度地提高轨控效率且避免系统角动量饱和现象.所提出的方法适用于喷气欠驱动控制系统，方法的有效性通过了数学仿真验证及在轨型号应用.     

An improved numerical method for constructing Halo/Lissajous orbits in a full solar system model

An improved numerical method that can construct Halo/Lissajous orbits in the vicinity of collinear libration points in a full solar system model is investigated. A full solar system gravitational model in the geocentric rotating coordinate system with a clear presentation of the angular velocity relative to the inertial coordinate system is proposed. An alternative way to determine patch points in the multiple shooting method is provided based on a dynamical analysis with Poincare′sections. By employing the new patch points and sequential quadratic programming, Halo orbits for L_1, L_2, and L_3 points as well as Lissajous orbits for L_1 and L_2 points in the EarthMoon system are generated with the proposed full solar system gravitational model to verify the effectiveness of the proposed method.